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A quantitative criterion to prove and analyze convergence within the numerical renormalization 
group (NRG) is introduced. By tracing out a few further NRG shells, the resulting reduced density 
matrices carry relevant information on numerical accuracy as well as entanglement. Their spectra 
can be analyzed twofold. The smallest eigenvalues provide a sensitive estimate of how much weight is 
discarded in the low energy description of latter iterations. As such, the discarded weight indicates 
in a site-specific manner whether sufficiently many states have been kept in a single NRG run. 
The largest eigenvalues of the reduced density matrices, on the other hand, lend themselves to a 
straightforward analysis in terms of entanglement spectra, which can be combined into entanglement 
flow diagrams. The latter show strong similarities with the well-known standard energy flow diagram 
of the NRG, supporting the prevalent usage of entanglement spectra to characterize different physical 
regimes. 
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I. INTRODUCTION 

The numerical renormalization group (NRG) 1 is a 
powerful method that provides a highly systematic non- 
perturbative approach to the wide realm of so-called 
quantum impurity systems. These consist of an ar- 
bitrary small quantum system (the impurity) in con- 
tact with a macroscopic non-interacting usually fermionic 
bath. Each part is simple to solve exactly on its 
own. In presence of interaction at the location of the 
impurity, however, the combination of both gives rise 
to strongly-correlated quantum-many-body phenomena." 
Wilson's logarithmic coarse-graining of the bath leads to 
a semi-infinite chain with exponentially decaying cou- 
plings, which justifies the concept of energy scale sep- 
aration. That is the Wilson chain can be diagonalized 
iteratively by adding one site at a time and retaining the 
lowest Mk states only. The obvious question, however, is 
how many states should one keep on average for conver- 
gence in this procedure? At a given iteration there is no 
quantitative a priori measure that indicates how many 
low-energy states are required for a proper description of 
the remaining low-energy physics. Usually, the only way 
to check convergence within the NRG is by repeating the 
entire calculation and showing that the results no longer 
change when further increasing Mk. Therefore an NRG 
calculation is typically run somewhat blindly for some 
pre-determined Mk- 

This somewhat uncontrolled truncation in the NRG 
is in stark contrast to the situation in the density ma- 
trix renormalization group (DMRG). 3-5 DMRG is based 
on a (strictly) variational principle, and as such has a 
clean well-defined truncation of the state space for part 
of the system through the discarded weight in its reduced 
density matrix. 4 In contrast to the less suggestive plain 
number Mk of states kept, the discarded weight repre- 
sents a reliable quantitative measure for the accuracy of 
a calculation. Within the DMRG, Mk can be easily ad- 



justed according to some predefined threshold in the dis- 
carded weight, instead. Motivated by DMRG then, an 
approximate similar criterion can be established within 
the NRG as will be shown in the following. The analysis 
requires a slightly longer chain, as shown schematically 
in Fig. 1. With the extra n$ sites traced out again from 
the ground state space of the enlarged system, this al- 
lows to estimate the discarded weight. The latter offers 
a quantitative convergence measure that is specifically of 
interest for numerically expensive models such as multi- 
channel models, or models where the energy scale sepa- 
ration along the Wilson chain might be in question due 
to modifications in the discretized Hamiltonian. In either 
case, a small discarded weight provides a strong indica- 
tion for converged NRG data. 

Furthermore, the reduced density matrices generated 
for the evaluation of the discarded weight also allow a 
quite different analysis in terms of their dominant cor- 
relations. In particular, combining their entanglement 
spectra into entanglement flow diagrams, this offers a 
complementary view to the usual NRG energy flow di- 
agram, which is entirely based on the analysis of the low- 
energy state space of a prior NRG run. 

The paper is thus organized as follows. In Sec. I the 
essentials of the numerical renormalization group are re- 




Figure 1: Schematic depiction of tracing out the low-energy 
sector of the Wilson chain at iteration n by including and 
analyzing no more NRG iterations. The impurity (dot) is 
entirely contained in the first site, while the bath is coarse- 
grained and mapped onto the remaining semi-infinite tight- 
binding chain of sites n = 0, 1, 2, 
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visited, including the construction of reduced density ma- 
trices. Sec. II then uses a specific set of reduced density 
matrices in the definition and analysis of the discarded 
weight within the NRG. Sec. Ill offers a complementary 
view on these reduced density matrices by analyzing their 
entanglement content in terms of entanglement spectra. 
Sec. IV, finally, summarizes and includes an outlook. 

A. The numerical renormalization group 

Within the NRG, the continuum of the non-interacting 
bath of half-bandwidth W is logarithmically coarse- 
grained in energy space, followed by an exact mapping 
onto a semi- infinite so-called Wilson-chain. 1 ' 6 The im- 
purity space is coupled to the first site of this chain 
only, as depicted schematically in Fig. 1. The logarith- 
mic coarse-graining is defined through the dimension- 
less discretization parameter A > 1. With the chem- 
ical potential at energy zero, the continuum of states 
in the energy intervals ±W[A~( n ~ z+1 \ A~(™~ z )] is ef- 
fectively represented by single fermionic levels (coarse- 
graining), including an arbitrary z-shift with z 6 [0, 1[. 7-9 
The subsequent exact mapping onto the semi-infinite 
chain (Lanczos tridiagonalization) 10 results in an effec- 
tive tight-binding chain with the exponentially decaying 
hopping t n ~ A - ™/ 2 between sites n and n + 1. For suf- 
ficiently large A, typically A > 1.7, this then justifies the 
essential NRG assumption of energy scale separation: by 
iterative diagonalization of the Wilson chain by adding 
one site at a time, large energies are considered first, 
with the (approximate) eigenstates at large energies dis- 
carded and considered unimportant in the description of 
the lower energy scales still to follow. Thus each site of 
the Wilson chain corresponds to an energy shell with a 
characteristic energy scale, 

^^^WA-l (1) 

Here, the prefactor was chosen such, that the rescaled 
couplings lim„_ i . 00 (t n /us n ) — 1 quickly approach unity 
for longer Wilson chains for arbitrary A and z-shift, with 
the discretization following the prescription of [9] for a 
flat hybridization, i.e. T(E) = T6(W - \E\). 

With H n the full Hamiltonian H of the Wilson chain 
up to and including site n, its low-energy eigenstates are 
given by the NRG eigenstates H n \s n ) = E™\s n ). Com- 
plemented by an arbitrary state \e n ) for the remainder 
of the system following site n, the NRG assumption of 
energy scale separation can be summarized then in the 
following approximation, 11 

H\se) n ~ E™\se) n , (2) 

i.e. the states \se) n = \s n ) <g> \e n ) are, to a good approxi- 
mation, also eigenstates of the entire Wilson chain. The 
energies E™ at iteration n are usually expressed relative 
to the ground state energy of that iteration, and rescaled 
by a factor ^ (A + 1) A - " 7 ' 2 cx u) n to resolve the energy 



shell at iteration n. The resulting energies are referred 
to as rescaled energies. For fully fermionic systems, they 
typically show an intrinsic even-odd behavior. Thus com- 
bining the rescaled energies vs. even and odd iterations 
n separately, this results in the standard energy flow di- 
agrams of the NRG. 1 ' 6 

The approximate many-body eigenstates |se) are con- 
structed iteratively, and therefore described in terms of 
matrix-product-states. 5,12-14 Each iterative step results 
in a basis transformation, encoded in an A-tensor, that 
combines an existing effective basis |s n ) for the system 
up to and including site n with the state space \o~) of site 
n + l, 

|s n +i) = V \s' n ,a n+ i)(s' n ,a n+1 \s n+1 ), (3) 

, S v ' 

with \s' n ,a n+ i) = \s' n ) (g> |u n+ i}. The orthogonality of 
state spaces, (s n +i\s' n+1 ) — 6 SS >, directly implies the or- 
thonormality relation for A-tensors, 4 

^K+dUK+d = 1. (4) 

Without truncation, the dimension M„ of the state space 
|s„) increases exponentially with the number of sites in- 
cluded, M n ~ d n , with d the dimension of a local Wilson 
site. Therefore the maximum number of states Mk, that 
one can maintain in a calculation, is quickly reached after 
n ~ log(A/K)/log(ci) iterations. For every subsequent 
iteration, the state space \s n ) is truncated by retaining 
the lowest Mk states in energy only. This leads to the 
distinction between \s%) and |s°) for kept and discarded 
states at iteration n, respectively. Correspondingly, this 
also splits the A-tensor into two parts, A^k and Akd-, 
that propagate the state kept space from the previous it- 
eration into the newly generated kept or discarded space, 
respectively. 

The truncation criteria with respect to a fixed pre- 
specified Mk can be softened towards a energy cutoff, 2 
i?K, that is taken constant in rescaled energies. For a 
fair comparison for different z-shifts, it will be specified 
in units of the energy scale w„ in Eq. (1). Since NRG 
data typically appears bunched at certain energies (e.g. 
see Fig. 3 later), Ek may hit a "gap" in the NRG spec- 
trum at some iteration, and the last "bunch" of states 
included may lie, on average, at clearly smaller energies 
than £k- Given the empirical importance of the first few 
NRG iterations, therefore as a safety measure, by default, 
£k was taken by 20% larger for the very first iteration 
where truncation occurred, i.e. using 1.2Ek there with 
Ek specified in context. Typical values are in the range 
E K = 5 . . . 8. 

The model system considered in this paper is the 
well-known standard single impurity Anderson model 
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(SIAM), 

h siam = e da fi a + Un dt h dl + £ ^/f (d\f 0a + /i.e.) 

+EE t »(tw+ k ) ( 5 ) 

(T n— 

with the operators dj. creating a particle with spin 

a G {f, -[•} at the impurity (at site n in the bath) , 
respectively, having hdc = a^a^. The energy Sda = 
Ed — ^-(fid-f — ftdi) is the spin dependent level-position of 
the impurity in the presence of a magnetic field B. Fur- 
thermore, U is the onsite Coulomb interaction and T the 
hybridization of the impurity with the bath. All parame- 
ters will be specified in units of the bandwidth W := 1 in 
context with the figure panels. The bath in Eq. (5) is al- 
ready represented in terms of a Wilson chain, 1 described 
by the semi-infinite tight binding chain (N — > oo) with 
exponentially decaying hopping amplitudes t„ ~ A~™/ 2 . 
In practice, N can be taken finite, with H n describing 
the Wilson chain up to and including site n < N. 

Charge and spin are conserved in the SIAM in Eq. (5), 
where, however, only the abelian part of the symmetries 
is included in the calculations. Hence the number of 
states Mk directly refers to the actual number of states 
kept in a calculation (in contrast to the dimension of 
reduced multiplet spaces with non-abelian symmetries). 
Similarly, also the discussion of the entanglement spec- 
tra further below will refer to the abelian symmetry la- 
bels which also applies when non-abelian symmetries are 
broken. Note that while, in general, a particle-hole sym- 
metric impurity setting will be used, this can be easily 
broken by applying a (small) gating potential to the im- 
purity level. Moreover, the SU(2) spin symmetry, in fact, 
will be broken explicitly by the application of an external 
magnetic field. 



B. Density matrices 

The NRG eigenbasis of Eq. (2) with respect to the dis- 
carded space forms a complete many-body eigenbasis. 11 
Initially introduced for the feat of real time-evolution 
within the NRG, this eigenbasis is actually applicable and 
tractable more generally within the NRG framework. 15 
In particular, this allows the clean calculation of correla- 
tion functions in terms of the full density matrix (FDM) 
in the many-body eigenbasis, 12 in that 



p(T) = U 



-m 



IE' 



-fSE" 



(6) 



with j3 = 1/kgT for arbitrary temperatures T, using 
non-rescaled energies E™ relative to a common energy 
reference, by construction of a thermal density matrix. 
Eq. (6) can be rewritten as p(T) = ^2 n w n (T)p n (T), i.e. 
a normalized distribution w n = 1 of the density ma- 
trices p n (T) generated in the basis of iteration n. 12 For 
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Figure 2: Backward update of given density matrix p n at 
iteration n. Blocks represent data spaces, lines correspond to 
indices. The lines connecting different blocks are contracted 
indices (i.e. indices summed over), such as a n , s n and s' n , 
while open lines represent open indices, e.g. the indices s n _i 
and si_i . 



a given temperature T, the distribution w n is strongly 
peaked around iteration nx that corresponds to the en- 
ergy scale of temperature. Hence temperature essentially 
terminates the Wilson chain. 

In this paper, however, mainly reduced density matri- 
ces derived from ground states will be considered, hence 
temperature is essentially zero. More generally then, con- 
sider an arbitrary density matrix defined in the many- 
body basis \s n ) of iteration n in either kept or discarded 
space, X e {K, D}, 



p[x] = 

rn — 



,[X] 



n s'£X 



(7) 



where p n (i.e. without the hat) represents the space of 
fxi 

matrix elements p s s , . The prototypical and well-known 
operation on such a density matrix is tracing out the last 
site n, 11 ' 12 ' 15 - 17 
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,KKt) 



(8) 



and 



written as a matrix product of the matrices A ' KX 

p n in the first line. Equation (8), in the following re- 
ferred to as backward update, introduces the notational 
shorthand V n for the bilinear product of the A- and A*- 
tensor at site n, that acts as a linear superoperator on the 
density matrix p n . The corresponding contraction pat- 
tern is shown in a simple graphical depiction in Fig. 2. 
By construction, the backward update of a density ma- 
trix in Eq. (8) always results in a density matrix in the 
kept space of the earlier iteration, and with Eq. (4) rep- 
resenting a complete positive map, Eq. (8) clearly also 
preserves the properties of a density matrix. 
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II. DISCARDED WEIGHT WITHIN THE NRG 

The standard notion of NRG is that it zooms in to- 
wards the low energy sector of a given many-body Hamil- 
tonian, while iteratively discarding states at higher ener- 
gies. Having a semi- infinite chain, this can continue to 
arbitrarily small energy scales, which enables NRG to re- 
solve dynamically generated small energy scales as they 
appear, for example, in the context of Kondo physics. 
From a variational point of view for matrix-product- 
states, this implies that the cost function can be iden- 
tified as 

lim (s N \H N \s N ) -> MIN, (9) 

N— >oo 

yielding the ground state |0)oo of the semi- infinite Wil- 
son chain. For a sufficiently long chain of total length 
N then included in a given calculation, the state |0)jv 
will be referred to as the overall ground state of this Wil- 
son chain. In fact, the cost function in Eq. (9) is well 
captured within the NRG through its principle of energy 
scale separation. 18 

If at a given iteration within the NRG states essen- 
tially decouple with respect to the low energy state space 
still to follow, these states will quickly and efficiently be 
discarded as high energy states. The truncation towards 
the low-energy sector also implies, that the state space at 
large energies is necessarily more crudely resolved, con- 
sistent with the coarser discretization there. The lowest 
Mk states kept at a given iteration n then are important 
for the correct description of the low-energy sector still 
to come. However, there is no real quantitative a-priori 
measure to indicate whether the number Mk of states to 
be kept is appropriate. Conversely, however, at a given 
iteration n one can ask whether all states kept a few it- 
erations earlier were actually important. This question 
can be answered entirely within the kept spaces of these 
iterations, hence is numerically cheap to analyze. 



A. Construction of reduced density matrices 

Consider the actual ground state space G at some arbi- 
trary but fixed iteration n' . In general, it may be g„<-fold 
degenerate, hence consider its fully mixed density matrix, 

po,^=mMM- ( io ) 

By construction, the number of eigenvalues of po tTl > un- 
equal zero, i.e. its Schmidt rank, is equal to g n i. Now, 
tracing out the last iteration nf, i.e. the lowest en- 
ergy scale included in po.n', is equivalent to the back- 
propagation /5q™ = V n 'Po, n ' in Eq. (8). Through this 
operation, the Schmidt rank will rise, in general, by a 
factor of d, with d the state space dimension of a Wil- 
son site. Repeating this process iteratively, this allows to 



trace out the no smallest energy shells in /3o, n ' • Thus with 
n' = n + no, this leads to the reduced density matrix, 

n+n 
l=n+l 

ss' 

which, by construction, is defined in the kept space of 
iteration n. The Schmidt rank will grow quickly, i.e. 
exponentially, in this process, until after n iterations, 
with 

n > ceil[log(M K )/ log(d)] (n « N), (12) 

it reaches the full dimension Mk of the kept space. Typ- 
ically, hq is much smaller compared to the full length N 
of the Wilson chain considered, and conversely also spec- 
ifies the initial number of NRG iterations in a forward 
direction that can be typically performed without trun- 
cation. For the definition of the discarded weight below, 
it is sufficient to stop the back-propagation of po,n+n a t 
this point. 

The reduced density matrix p n,n °l generated in 
Eq. (11) is, in general, not diagonal in the energy eigen- 
basis |s„), since through the traced out lower-energy sites 
it does know about an enlarged system. Its eigenvectors 
are described by a unitary transformation iti™'™ ' within 
the NRG eigenstates kept at iteration n, 

Kn )^£4r ol K K >, 

s' 

With P t' n \n-,n )=p [ r l0] \rn,n ), (13) 

where the index r shall refer to the eigenstates of the 
reduced density matrix, in contrast to the index s for the 
energy eigenstates. Here, the eigenvalue pi™'™ ' describes 
the importance of a specific linear superposition of NRG 
eigenstates at iteration n for the low-energy description 
of latter iterations. 

This offers two routes for the analysis of the density 
matrices /5 ™'™ '. (i) Adhering to the energy eigenbasis 
of the NRG, the importance of the kept state \s%) at 
eigenenergy E™ for the latter low-energy physics is given 
by the expectation value 

pl^l^l^" 01 !^}, (14) 

i.e. the diagonal matrix elements pi™'™ '. Alternatively, 
(ii) using the eigenbasis of the reduced density matrices, 
the weights of these states are given by the eigenvalues 
pi™'™ ', while now their energies are given by the expec- 
tation values 

E^ = {r n , no \H n \r n , no ). (15) 

Both routes will be analyzed and compared in the fol- 
lowing. However, the actual eigendecomposition of the 
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Figure 3: (Color online) Weight distribution of energy eigen- 
states over full NRG run at fixed M K = 512 for the SIAM 
[Eq. (5): U = 0.20, e d = -U/2, V = 0.01]. The main panel 
shows the rescaled eigenenergies E™ vs. their weights p!,™'™ ' 
as in Eq. (14). Data is shown only for those iterations where 
truncation occurred, with data from the same iteration shown 
in the same color. The two iterations with smallest (largest) 
energy range, n max (n m i n ), are highlighted in strong colors 
(black diamonds (red crosses)), respectively, while light col- 
ors are used for all other iterations. The top [right] panel 
shows the energy [weight] distribution u(E) [f{p)], Eq. (16) 
[Eq. (17)], respectively, for the data in the main panel, with 
matching energy E [weight p] axis. The binning referred to 
in the text to Eqs. (16) and (17) is indicated by the intervals 
between the data points in the top and right panel. 



(ii) the average distribution of the weights p^ 

N , 

Kp) = ^E E i. ( 17 ) 

n=l p<pi<p+dp 

and (iii) their average dependence on each other 

N , 

p( e ) = wseJ2 E Pi ( 18a ) 

n E<Ei<E+dE 

£ Ke~ KE . (18b) 

Here some appropriate linear (logarithmic) binning of 
the data is assumed with energy (weight) intervals dE 
(dp), respectively. In particular, the densities in Eqs. (16) 
and (17) are clearly dependent on these binning inter- 
vals, which therefore will properly indicated in the sub- 
sequent plots. The prime in the summation and the nor- 
malization indicates that only those iterations n are in- 
cluded where state space truncation occurred, i.e. typ- 
ically n > uq. The total number of these iterations is 
given by N'. With chosen normalization then, the sum 
over the binned v(E) and v(p) data both yield the av- 
erage number of kept states, while the integrated weight 
distribution p(E) in Eq. (18a) is normalized to 1, since 
tr(p) ~ J °° p{E) dE = 1. As will be seen later, the 
weight distribution p(E) typically shows a clear exponen- 
tial decay with a characteristic exponent K , cLS indicated 
already in Eq. (18b), with the prefactor chosen such that 
it also preserves normalization. 



1. Energy eigenbasis 



reduced density matrices will be preferred for the remain- 
der of the paper as explained. 

In either case, a set of states i with (average) energy Ei 
is given together with their respective (average) weight p± 
that represents the states importance for latter iterations. 
For the first [second] route above this data is given by 

(E^,p [ r io] ) [(£j" ; ™ o] ,p[ n;nQl )], respectively. Given that 
the reduced density matrix /5q"'™ ', by construction, ex- 
ists in the kept space only, therefore also all states i re- 
fer to the kept space or a linear superpositions thereof. 
Moreover, for every iteration, the weights pi are normal- 
ized, i.e. they are positive and add up to 1, while by 
combining data from different iterations, the energies Ei 
are always specified in rescaled units. 

The resulting data (Ei,pi) then is clearly correlated. 
It is analyzed threefold, (i) in terms of the average dis- 
tribution of the rescaled energies 2?, 



K£) = ^E' E l ( 16 ) 

n=l E<Ei<E+dE 



The correlation between the eigenenergies £J™ and their 
corresponding weights pi"'™ ' is plotted as a scatter plot 
in the main panel of Fig. 3. The model analyzed is the 
SIAM in Eq. (5) in the Kondo regime using a fixed num- 
ber of kept states, with all parameters specified in the 
figure caption. The weights pi™'™ ' clearly diminish expo- 
nentially with energy, which is intuitively expected as a 
consequence of energy scale separation within the NRG. 
The integrated weight distribution p(E) (dashed black 
line, cf. Eq. 18a), shows a clear exponential decay with 
an exponent k ~ 2.7. As seen in Fig. 3, this distribu- 
tion clearly also serves as an upper bound of the weights 

[n;no] , 

p s at a given energy. 

The upper panel in Fig. 3 shows the distribution v(E) 
in Eq. (16) of the energies E™ plotted in the main panel 
(matching horizontal axis). This distribution shows a 
strong increase with energy E, consistent with the no- 
tion that the many-body phase space grows quickly as 
the available energy for excitations becomes larger. To- 
wards large energies, eventually, the data is necessarily 
truncated to the finite number Mk of kept states, which 
leads to a drop in the density v(E). The exact bound- 
ary with respect to energy is somewhat blurred, though, 
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since in given case fixed Mk allows the energy range to 
vary for different iterations n. The right panel of Fig. 3, 
on the other hand, shows the distribution v{p) in Eq. (17) 

of the weights pj™'™ ' plotted in the main panel (match- 
ing vertical axis) . This distribution is peaked around the 
largest weights pi"'™ ' for the largest energies E™. 

The data in the main panel of Fig. 3 is typically 
bunched around a set of energies for a fixed iteration 
n. This is also reflected in the distribution v (E) in the 
upper panel of Fig. 3, and is due to the discretization of 
the model. Moreover, two iterations are highlighted in 
strong colors. These correspond to the iterations whose 
energy range is smallest (n min = 6, red bullets) or largest 
("■max = 74, black diamonds). Intuitively, the largest nu- 
merical error is expected from iterations such as n m i n (red 
bullets) since through Eq. (18b), stopping at premature 
energies directly translates to largest missing, i.e. dis- 
carded weight in the density matrix. As an aside, this 
serves as a strong argument in favor of truncation w.r.t. 
a fixed energy cutoff rather than a fixed number Mk 
of states. Fixed Ek, however, also introduces more noise 
to the data in particular for higher lying states. Hence 
both truncations will be used and pointed out in context. 

The weights pi™'™ ' in the main panel of Fig. 3 show 
significant vertical spread, which translates into a pro- 
nounced tail towards exponentially smaller p in the dis- 
tribution v{p) in the right panel. For a given energy E 
therefore, many of the states have order of magnitudes 
lower weight than the top-most weights close to p(E) in 
the main panel. This indicates that the energy repre- 
sentation with its corresponding diagonal weights pi™'™ ' 
is not necessarily the optimal basis to analyze accuracy. 
Moreover, note that using the energy eigenbasis \s n ) with 
energies E™ in the analysis of the reduced density ma- 
trices, this actually mingles the energy scales of an ef- 
fectively larger system H n+no with the basis generated 
w.r.t. H n only. 



2. Eigenbasis of reduced density matrices 

From the point of view of a variationally optimal repre- 
sentation of the ground state space of an enlarged system, 
on the other hand, one is directly led to the eigenspec- 
trum of the reduced density matrix, as exemplified within 
DMRG. 3 The analysis of Fig. 3 therefore is repeated for 
the same underlying Wilson chain, yet with two modifica- 
tions: (i) the eigendecomposition of the reduced density 
matrices in Eq. (13) together with Eq. (15) is used in- 
stead of the energy eigenbasis, and furthermore (ii) the 
NRG truncation criterion is based on a fixed energy cut- 
off, -Ek = 6. The results are shown in Fig. 4, with striking 
quantitative differences compared to Fig. 3. The spread 
in the scatter plot is significantly narrowed, and overall, 
the data decays much faster with k ~ 4.6, cf. Eq. (18b). 
Therefore this leads to a clearly improved separation of 
the actually relevant states for the subsequent description 




Figure 4: (Color online) Similar analysis as in Fig. 3 (see 
caption there for further information) for the same underlying 
Hamiltonian, except that the eigenspectrum of the reduced 
density matrices in Eq. (13) was used together with Eq. (15) 
and a fixed energy cutoff Ek ~ 6. Similar to Fig. 3, only those 
iterations are shown where truncation occurred (same color 
for data from the same iteration), with the same two iterations 
highlighted as in Fig. 3, indicated by n± and ri2. The estimate 
for the overall discarded weig ht e° =5% =i 6 • 10" 12 as defined 
in Eq. (21) is indicated by the horizontal dashed line. 



of the lower-energy scales. This suggests that many of 
the NRG eigenstates, as their energy increases, loose im- 
portance much faster as compared to Fig. 3, despite the 
relatively large diagonal weights p s in the density matrix 
still seen there. In a sense, the weights there represent 
mere matrix-elements in a non-diagonal representation. 

The iterations highlighted in Fig. 4 are the same iter- 
ations as in Fig. 3. Given a fixed energy cutoff E-& = 6 
here, however, both have a comparable energy range 
(hence the altered notation ni and n%), with the num- 
ber Mk of kept states varying from ~ 1000 at very early 
iterations (in particular iteration ni), down to ~ 250 
at late iterations (such as iteration 712). Note also the 
markedly fewer data points seen for iteration 712 . This is 
only partly due to the reduced number of states, as there 
are also large systematic (approximate) degeneracies at 
the strong-coupling Kondo fixed point already reached 
at this iteration. This results in many of the black dia- 
monds lying indistinguishably on top of each other (see 
also discussion on entanglement spectra later). 

As seen from above discussion, rather than taking the 
energy eigenstates |s n ) and the corresponding diagonal 
matrix elements pi™'™ ' (Fig. 3), the eigenvalues p[™'™ ' of 

the reduced density matrix />q™'™ ' do represent a clearly 
better choice for the analysis of accuracy or entanglement 
in the system (Fig. 4), and thus will be used henceforth. 
This prescription also shows a more systematic exponen- 
tial decay all the way down to numerical double precision 
noise (10 -16 ), with the decay rate k of p(E) roughly in- 
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Figure 5: (Color online) Discarded weig ht e D for the SIAM 
[Eq. (5): U = 0.20, e d = -(7/2, F = 0.01 (same parame- 
ters as in Fig. 3), with T K oi 1.23 ■ 10~ 5 ]. Panel (a) shows 
the discarded weight ef x ) defined in Eq. (21) vs. My. us- 
ing no £ {6,7,8}. The data £°=5% is shown in solid lines, 
while the data based on the minimum eigenvalue of p^'™ ' 
(cf. Eq. 19b) is shown in dashed lines. The distribution of 
the discarded weight e„ x along the Wilson chain is shown in 
the inset for M K G {128,256,512, 1024}, also marked by the 
vertical dashed lines in the main panels. Panel (b) shows the 
conductance go vs. Mk in units of 2e 2 /h while using a set of 
shifted discretizations, with the z-values as specified. Conver- 
gence in the conductance towards the expected unitary limit 
is seen for M K > 400, i.e. e° < 10" 12 . 



dependent of the discretization parameter A. 



3. Definition of discarded weight 

With the motivation above, the definition of the dis- 
carded weight is based on the eigendecomposition of the 
reduced density matrices p n ' n °^ in Eq. (11), using the 
combined data of Eq. (13) and Eq. (15). In terms of 
Fig. 4, adding more states to the calculation essentially 
extends the data to larger energies and smaller weights, 
while the large-weight low-energy sector already remains 
widely intact. Therefore the largest discarded weight, i.e. 
the weight missing by states not included and hence not 
available, can be estimated, to a good approximation, up 
to an overall prefactor by the smallest weights in the kept 
state space, which are easily accessible. Given the expo- 



nential decay of the weights together with the residual 
spread in the data as seen in Fig. 4, the discarded weight 
at given iteration n can thus be defined through the av- 
erage weights pr l,n °^ for the highest energies i^j™'™ ' in 
the kept space, 

£» x nn = (/4" ; " o] / , , ,„„• (19a) 

The parameter x <^ 1 is considered small, yet is cho- 
sen large enough (typically x — 0.05) to average over the 
residual spread of weights. Alternatively and for compar- 
ison, an even simpler measure in terms of the minimum 
eigenvalue of Pq™'™ ' will be considered, 

4„ ^min(p[" : "«]), (19b) 

which no longer makes any explicit reference to energies. 
Note that even though e^n or e^. n , written £ m; ( ^ in 
short, are purely determined within the kept space, they 
clearly represent a sensible estimate for the discarded 
weight at iteration n, i.e. e n (x) ~ e„'4, defined as the 
fraction of relevant state space missing from the latter 
description of the low energy physics. If no truncation 
has occurred at iteration n, however, such as typically 
for the first n < uq iterations, of course, then there is no 
truncation error either, hence £ n {x) = for these itera- 
tions. 

In summary, the discarded weight at iteration n is 
defined as follows, 

£ D (x) — f Enjno m tne presence of truncation 

[0 without truncation at iteration n. 

(20) 

Here e n -nl can be determined efficiently by including 

and analyzing uq further NRG iterations within the kept 

space, where typically no <C N, cf. Eq. (12). The overall 

discarded weight eP s of a full NRG run then is taken, for 
° (x) _ 

simplicity, as the largest discarded weight per iteration, 
ef x) ^max( £ >>), (21) 

Using x = 5% as in Eq. (19a), the discarded weight for 
the NRG run in Fig. 4 is estimated by ~ 6-1CT 12 , indi- 
cated by the horizontal dashed line. As seen from Fig. 4, 
the overall discarded weig ht e° for an NRG run essen- 
tially coincides with p(E) at the largest energies within 
the kept space. On the other hand, e D , i.e. without the 
usage of x based on the plain minimum eigenvalue of the 
reduced density matrices p Q n ' n °\ cf. Eq. (19b), will in 
general lie a (constant) few orders of magnitude lower, 
as it happens, for example, for the data in Fig. 4. Never- 
theless, as will be shown in the following, up to an overall 
global prefactor the discarded weight based on either, e D 
or eP, both behave in an essentially similar fashion. 



B. Application 

The discarded weight ePj defined in Eq. (21) sensi- 
tively depends on the number Mk of states kept or the 
energy threshold Ek- From Fig. 4 one expects a strongly 
diminishing discarded weight with increasing Mk or Ek, 
a quantitative analysis of which is presented in Figs. 5 
and 6 for the SIAM. Figure 5 analyzes the dependence 
of the discarded weight e? * on the number Mk of states 

kept. As seen in panel (a), the discarded weight ePx 
strongly decays with Mk, with minor variations when 
a new Wilson shell is fully included without truncation, 
e.g. at M K € {256, 1024}. With panel (a) being a log-log 
plot, the decay of the discarded weight with Mk rather 
resembles a polynomial convergence, yet with very large 
power (on the order of 10). The reason for the slower 
than exponential decay is due to the strong increase in 
the density of states v{E) of the full many-body eigen- 
spectrum with increasing E as discussed with Figs. 3 and 
4. 

Together with the analysis of the discarded weight in 
Fig. 5, an independent physical check for convergence is 
provided by the numerically computed conductance go 
in units of 2e 2 /h shown in Fig. 5(b). The conductance 
was calculated via the (spin-resolved) spectral function 
A (tr) (w) = / §e i " t {{d a {t),dt}) T of the impurity level, 
with g = irT J dcj(— I^)A(cj). Here the Fermi func- 
tion f(uj) and the spectral function A(w) are evaluated at 
small but finite temperature T ~ 6 • 10~ 8 , which is much 
smaller than the Kondo temperature of 7k — 1.23 • 10~ 5 
for given parameter set and corresponds to the energy 
scale close to the end of the Wilson chain, having A = 2 
and N = 60. Expecting go = 1 for the symmetric SIAM, 
the data in Fig. 5 indicates convergence for Mk > 400. 
The data for smaller Mk is not yet converged, and there- 
fore (strongly) depends on numerical details, such as non- 
averaged z-shifts. 8 ' 9 

With Mk being constant, the energy of the topmost 
kept states can vary significantly with Wilson shell n, 
which directly also leads to a clear dependence of the 
discarded weight eP * on n. This is shown in the inset 
to panel (a) for the set of different values of Mk marked 
in the main panels by the vertical dashed lines. The dis- 
carded weight clearly varies over more than three or- 
ders of magnitude within a single NRG run, irrespective 
of the actual Mk. In particular, one can see that earlier 
iterations dominate the discarded weight e° for physical 
reasons. In the strong-coupling regime for n > tik (with 
iteration % — 35 corresponding to the energy scale of 
Tk), the discarded weight is smallest, while for the inter- 
mediate free orbital or local moment regime for n < tik, 
these regimes require a larger number of states for com- 
parable numerical accuracy from a physical point of view, 
indeed. 

Given the underlying energy scale separation of the 
NRG, a straightforward way to obtain a more equally 
distributed e n <x) is achieved using an energy cutoff Ek, 




Figure 6: (Color online) Similar analysis as in Fig. 5, yet for 
truncation with respect to fixed energy Ek- For several values 
of Ek marked by the vertical dashed lines in the main panels, 
the distribution of the discarded weight £^ x along the Wilson 
shell n is shown in the inset to panel (a). With Mk allowed 
to vary over a wider range, panel (c) shows the correlation of 
Mk with _Ek, plotting average, minimum, and maximum of 
Mk along the Wilson chain. For the average Mk, data for 
different z-shifts is shown (several lines on top of each other, 
with same color coding as in panel b). 

as demonstrated in Fig. 6 for exactly the same system 
as in Fig. 5 otherwise. For the values of Ek indicated 
by the vertical dashed lines in the main panels, the in- 
set to Fig. 6(a) shows the distribution of En x ■ By con- 
struction, the discarded weight is, up to even-odd oscilla- 
tions, clearly more uniformly distributed over the Wilson 
shells as compared to the case of fixed Mk in Fig. 5(a). 
The discarded weight in panel (a) clearly diminishes ex- 
ponentially with £?Ki y e t with pronounced intermedi- 
ate plateaus since the discrete eigenenergies within an 
NRG run are usually bunched around certain energies. 
The corresponding average Mk as function of Ek, never- 
theless, follows a rather smooth monotonic behavior, as 
shown in panel (c). Given fixed Ek, however, clear vari- 
ations of Mk are seen within a given NRG run, hence 
also smallest and largest Mk are shown in panel (c). Ig- 
noring iterations without truncation, in given example, 
typically the largest Mk is required at early iterations, 
while the smallest Mk are encountered in the strong cou- 
pling regime at late iterations n > tik- 

The calculated conductance shown in panel (b) con- 
verges clearly more uniformly with increasing Ek as 
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compared to Fig. 5(b). In particular, it indicates con- 
verged NRG data for Ek > 5.5, which corresponds to 
£ x ~ Therefore in both settings, for constant 

Mk in Fig. 5 as well as for constant Ek in Fig. 6, conver- 
gence of the physical data is found for a similar discarded 
weight of < 10~ 12 with a negligible dependence on uq. 
This value therefore is considered a sufficient bound in ac- 
curacy to capture the main physics, with other quantities 
such as the NRG energy flow diagram already also well 
converged. 

Alternatively, using the plain minimum of the eigen- 
values of the reduced density matrices in Eq. (19b), this 
leads to convergence for e D < 1CP 16 . Given that e D refers 
to the minimum eigenvalue in the kept space, e D consis- 
tently lies about three orders of magnitudes lower than 
and is considered a lower bound to the actual discarded 
weight. While e D fluctuates slightly more strongly com- 
pared to e° owing to the fact that it is not an averaged 
quantity such as e®, it nevertheless follows a similar con- 
sistent picture in terms of convergence with the number 
Mk of states kept or the energy Ek used for truncation. 
In this sense, either discarded weight, e D as well as e°, 
can be used quite generally as a quantitative measure, 
indeed, to demonstrate accuracy within the NRG. In or- 
der to avoid confusion, however, it shall be made clear 
which one is used. 



III. ENTANGLEMENT SPECTRA 

The reduced density matrices pg 1 '" ' clearly also carry 
physical information in terms of entanglement along the 
Wilson chain. This is provided by the high end of their 
spectral decomposition. There the exact details of the 
largest eigenvalues of /5 n '™ ' are of interest, which do vary 
with no over a wider range depending on the underlying 
physics. Hence, in the following, the actual entanglement 
spectra will be calculated with respect to the reduced 
density matrices p Q of the overall ground state of the 
system, 



Po 



= lim 

no— >-oc 



Pb 



n;n ] 



7i;N-n] 



(22) 



The length N of the Wilson chain is taken sufficiently 
large, such that the energy scale of the last iteration N 
is much smaller than any other energy scale in the sys- 
tem. Temperature is therefore essentially zero. For com- 
parison, also the truncated entanglement spectra will be 
calculated from /5 n ' n °' for finite small n n , with n spec- 
ified in context. Motivated by the discussion following 
Eq. (6), the latter analysis can be linked to finite tem- 
perature settings. 



1. General definition 

The partitioning of the Wilson chain into two parts, the 
chain up to and including site n (part A) , and the traced 
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Figure 7: (Color online) Comparison of the standard NRG 
energy flow diagram (left panels) to the entanglement flow 
diagram (right panels) for the symmetric SIAM [U = 0.2, 
E d = -(7/2, T = 0.01, T K = 1.2 ■ 10" 5 ; A = 2, M K = 512, 
N — 80], with top (bottom) panels for even (odd) iterations, 
respectively. In addition to the actual entanglement flow di- 
agram obtained from the ground state of the last iteration at 
TV = 80 (black lines), also the truncated entanglement flow di- 
agram is shown, using no — 8 (orange (gray) lines). For better 
comparison with the energy flow diagram, the entanglement 
spectra (right panels) are also shifted at every iteration with 
respect to the smallest entanglement energy min(£). The y- 
scale of the entanglement spectra was adjusted to best match 
the energy fixed point spectrum in the left panels. Degenera- 
cies of energies at large n, i.e. lines lying indistinguishably 
on top of each other, are specified by the numbers on top of 
the lines in all panels. 
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Figure 8: (Color online) Comparison of spin-resolved fixed 
point spectra for the symmetric SIAM in Fig. 7 in the SC 
regime (n = 60) . Panel (a) [b] show the energy [entanglement] 
fixed point spectrum, respectively, vs. spin symmetry quan- 
tum number S z . For all low-energy multiplets the underly- 
ing (approximate) degeneracy is indicated. The entanglement 
spectrum is shifted w.r.t. to its lowest energy and scaled to 
mAeTfTO3fl§F g offtop&te m ^ In par- 
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ticular, this allows to make use of the recently introduced 
entanglement spectra (ES) 19 for the physical characteri- 
zation of a given wave function. Here these entanglement 
spectra provide a powerful tool for the systematic anal- 
ysis of the physical correlations in the reduced density 
matrices jffl in Eq. (22). 

Consider a given wave function of a some system par- 
titioned into parts A and B. The reduced density matrix 
Pa = trs(/o) is obtained by tracing out part B of the 
overall density matrix p. Within this setting, the en- 
tanglement spectrum is defined as the spectrum of the 
fictitious Hamiltonian H^, 19 

p A =: exp(-Hf). 

One may assume an effective inverse temperature ft := 1 
in order to make contact with a thermal density matrix. 
This f3 also sets the (otherwise arbitrary) energy scale 
in the per se dimensionless H^. With pa a positive op- 
erator, the entanglement spectrum £ r is defined as the 
eigenvalues of H^, i.e. 

£ r :=-logp r , (23) 

with p r the spectral decomposition of the reduced den- 
sity matrix pa- Particular information can be read off 
from the entanglement spectrum as soon as there is a 
rich amount of quantum numbers specifying the entan- 
glement levels and when entanglement gaps appear which 
separate a low-lying generic set of levels from irrelevant 
background correlations. 19-21 The spectra p r and £ r are 
independent of whether A or B is traced out, while of 
course, they are dependent on the specific choice of the 
partitioning. For entanglement spectra, the partition- 
ing typically occurs in real space for gapped systems, 
analyzing the edge of the thus created boundary, while 
for gapless systems momentum space is preferred. 20 The 
second case then is consistent with the systematic NRG 
prescription of energy scales based on the underlying dis- 
cretization in energy (momentum) space. 

By construction, the dominant correlations between 
systems A and B correspond to the lowest entanglement 
energies £ r , while weaker correlations will rise to higher 
energies. By tracing out a major part of the system, en- 
tanglement spectra provide significantly more informa- 
tion, say, than just the entanglement entropy between A 
and B. In particular, it has been shown that it provides 
finger prints of the underlying physics, and as such al- 
lows to characterize the physical nature of a given wave 
function. 19,20 This analysis is therefore entirely targeted 
at a given (ground state) wave function, without any fur- 
ther reference to an underlying physical Hamiltonian that 
it may have originated from. 

2. Application to NRG 

The general concept of the entanglement spectra can 
be readily transferred to the NRG. At each iteration n, 



the reduced density matrix p [ Q in Eq. (22) is computed 
and diagonalized, with its eigenspectrum mapped onto 
the entanglement spectrum in Eq. (23). Collecting these 
spectra and plotting them vs. iteration index n for even 
and odd iterations separately, the result will be referred 
to as entanglement flow diagram, in complete analogy 
to the standard energy flow diagrams of the NRG. For 
comparison, also the truncated entanglement spectra for 
finite small no will be analyzed, which in their combina- 
tion will be referred to as truncated entanglement flow 
diagram. In either case, the entanglement spectra are 
obtained in a backward sweep, purely based on the iter- 
ative low-energy Hilbert-space decomposition of a prior 
NRG run in terms of the A-tensors in Eq. (3). This is in 
contrast to the energy flow diagram, which is calculated 
with increasing shell index n in a forward sweep making 
explicit reference to the Hamiltonian. 

The entanglement spectra were calculated for the sym- 
metric SIAM in the absence of magnetic field. The re- 
sulting entanglement flow diagram is presented in Fig. 7 
together with a direct comparison to the standard NRG 
energy flow diagram. The data is plotted for even (odd) 
Wilson shells n in the upper (lower) panels, respectively. 
The energy flow diagram, shown in the left panels, clearly 
distinguish the well-known physical regimes of the SIAM, 
namely the free orbital regime (FO; n < 10), the local 
moment regime (LM; 10 < n < «k), and the strong 
coupling regime (SC; n > ^k)> where % ~ 35 corre- 
sponds to the energy scale of the Kondo temperature 
T K = 1.2 • 10~ 5 , having A = 2. All degeneracies for 
n > ?ik are explicitly specified in Fig. 7. In particular, 
for even iterations, the ground state is unique through- 
out, e.g. the Kondo singlet for n > jik (panel a), while for 
odd-iterations the ground state space at small energies is 
four-fold degenerate due to the particle-hole symmetric 
parameter set (panel c). 

Interestingly, a very similar picture emerges from the 
entanglement flow diagram in right panels, Fig. 7(b+d) 
(black lines). For comparison, also the truncated entan- 
glement spectra are shown using n$ — 8 (orange (gray) 
lines), which in given case converge rapidly, in fact ex- 
ponentially, with increasing no < 10 towards the actual 
entanglement flow diagram. The latter then mimic the 
energy flow diagram in the left panels over a wide range. 
For example, the convergence towards the Kondo fixed 
point occurs around similar iterations, and even the de- 
generacies of the lowest states of the energy flow dia- 
gram are exactly recovered by the entanglement spectra. 
The latter also holds on the symmetry-resolved level, as 
demonstrated in Fig. 8 for the even iteration n — 60 (see 
later discussion). Nevertheless, looking more closely, a 
few notable qualitative differences of the entanglement 
flow diagrams in right panels of Fig. 7 are seen compared 
to the energy flow diagrams in the left panels. Overall, 
the entanglement flow diagrams appear shifted by about 
5 iterations to larger energies. This can be understood, 
considering that the entanglement spectra are calculated 
for enlarged systems together with the rapid convergence 
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Figure 9: (Color online) Comparison of the standard NRG 
energy flow diagram (left panels) to the entanglement flow 
diagram (right panels) for the SIAM at finite magnetic field 
(same analysis as in Fig. 7, otherwise, see caption there for 
details, with same model parameters, except B = 2 ■ 1Q~° ~ 
1.6 T K ). 
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Figure 10: (Color online) Comparison of spin-resolved fixed 
point spectra for the SIAM at finite magnetic field in Fig. 9 
at the even iteration n = 60 (similar analysis as in Fig. 8 
otherwise) . 



with increasing no in given case. However, there are fur- 
ther pronounced differences with the energy flow diagram 
for the earliest iterations in the FO regime, n < 10. 

These latter differences in the entanglement flow dia- 
gram can be significantly enhanced by turning on a mag- 
netic field on the order of the Kondo temperature, as 
shown in Fig. 9 for B = 1.6 Tk- This corresponds to the 
energy scale at iteration iib ~ 32, given A = 2. The 
magnetic field has been chosen such, that for late iter- 
ations n 3> jib the fixed point spectrum for even and 
odd iterations become essentially the same (compare the 



low-energy fixed point spectra in panels a(b) to c(d), re- 
spectively). Due to the magnetic field, the Kondo sin- 
glet (previously the unique state at even iterations) is 
largely destroyed for n > with a sizeable magneti- 
zation at the impurity. Clearly, the NRG eigenbasis at 
early iterations n < tib does not yet know about the 
small energy physics to come {e.g. the small B ~ Tj< ap- 
plied in given case). Therefore the energy flow diagram 
essentially remains unaltered there, when compared to 
the case without magnetic field in Fig. 7(a+c). The flow 
changes strongly only starting from the energy scale of 
the magnetic field value, i.e. for n > ijb where it moves 
into a different fixed-point spectrum. In particular, there 
also emerges a unique state now in the energy flow dia- 
gram for odd iterations for n ^> i.e. the symmetry 
broken spinful state favored by the magnetic field. By 
including magnetic field, the entanglement flow diagram 
shows pronounced differences from the energy flow dia- 
gram for n < 71b , which includes large portions of the LM 
regime. While the energy spectrum up to and including 
site n is ignorant of the low-energy physics to come, this 
very low-energy physics is captured by the reduced den- 
sity matrices and thus reflected in the entanglement flow 
diagram. 24 

Consider the entanglement spectra derived from the 
overall ground state (black lines) in Fig. 9(b+d). In 
panel (d) the ground state remains unique throughout, 
i.e. remembers the symmetry broken magnetic state, de- 
termined at much lower energy scales, all the way up to 
the largest energies. Within the split-up lowest energy 
space with subsequent degeneracies [1-2-1] in panel (d) 
for n ^ Tie (to be called [1-2-1] configuration), the first 
and second excited states cross each other with decreas- 
ing n leading to a [1-1-2] configuration for small n, i.e. 
large energies. Nevertheless, the singly degenerate ex- 
cited state clearly remains split-off, and does not merge 
with the ground state, which is in strong contrast to 
the energy flow diagram in panel (c) with a [2-2] con- 
figuration for n <C riB- This degeneracy in the ground 
state space that is ignorant of the small magnetic field is 
partly reflected only in the truncated entanglement flow 
diagram. Using small uq (orange (gray) lines in panels 
d), this eventually also misses the low energy physics. 
Therefore these spectra in panel (d) eventually are also 
in a [2-2] configuration for the smallest n, with a more 
irregular transient behavior with increasing n. A similar 
trend is also observed for even iterations in panels (a+b) . 
While the ground state remains unique for all iterations 
in both panels, the entanglement flow in panel (b) tends 
to split off the excited levels right above the lowest [1-2-1] 
state space configuration for small n. For the truncated 
entanglement flow, on the other hand, the lines of these 
excited levels remain entangled with higher excitations, 
which is similar to the situation in the energy flow dia- 
gram in panel (a). 

Nevertheless, the low-energy fixed-point spectra for 
n 7Jb again agree well for both the energy and entan- 
glement flow diagram in Fig. 9, which again also holds 
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for the symmetry-resolved spectra, as demonstrated for 
the even iteration n — 60 in Fig. 10. This agreement in 
the spectra of the stable low-energy fixed point, present 
in both the non-magnetic as well as the magnetic case, 
is understood as a generic feature. There both, the en- 
ergy eigenstates as well as the reduced density matrices 
are deeply rooted in the low-energy physics, i.e. of the 
overall ground state of the system at T —¥ 0, and hence 
present a consistent description of the system. 

The detailed structure of the energy fixed point spec- 
tra provides clear physical information. 1,2 This includes, 
for example, phase shifts if a Fermi-liquid point of view 
is supported as is the case for the SIAM. This then di- 
rectly explains all of the splittings and degeneracies in 
the low energy sector of the energy fixed point spectra. 
For example, consider the energy spectrum in Fig. 8(a) 
for the fully symmetric SIAM in the non-magnetic case. 
Note that while spin-resolved spectra are shown in Fig. 8, 
in given case the charge-resolved spectra would look ex- 
actly the same due to particle-hole symmetry. With the 
spectra shown for an even iteration, the ground state is 
unique, i.e. represents the Kondo singlet with S z = 0. 
The first excited states for S z — + \ , correspond to an ex- 
tra particle with spin-up or a hole with spin-down. Given 
particle-hole symmetry, both processes have the same en- 
ergy 8/2 = 0.63 (in rescaled energy units), and hence are 
two-fold degenerate, indicated by the number on top of 
the level in Fig. 8. By symmetry, the same excitations 
exist for 2S Z — — 1, leading to the [2-2] degeneracy (4 
states) in the lowest excitations in Fig. 8(a). The next 
higher excitation combines two of above processes. This 
leads to a total of 6 excitations, all with energy 8 and 
distributed over 2S Z E {—2,0, +2}. Here two of the ex- 
citations at 2S Z = correspond to the extraction or an- 
nihilation of two particles with opposite spin. This fully 
explains the [1-4-1] degeneracy of the excited states at 
energy 8 = 1.26 in Fig. 8(a), and also the combined 6- 
fold degeneracy seen in the energy flow diagram seen at 
this energy in Fig. 7(a) . The argument can be continued 
along similar lines to explain the [4-4] (8 states) and [4- 
9-4] (17 states) degenerate subspaces of the next higher 
excitations. Excitations with even higher energy eventu- 
ally have missing levels due to NRG truncation. 

The same analysis as for the energy spectra, however, 
cannot be applied with equal rigor to the entanglement 
spectra. While the ground state [1] and the lowest [2-2] 
and [1-4-1] excitations in Fig. 8(b) fully agree in sym- 
metries, degeneracy and also in the precise relative level 
spacing, the next higher [4-4] excitation in panel (a) is 
broken up in Fig. 8(b), with some of the levels shifting to 
higher entanglement energy. Nevertheless, the degener- 
ate set [2-10-10-2] further up in energy still again equally 
appears for both, energy and entanglement spectra. 

The same analysis as in Fig. 8, is repeated for the mag- 
netic case in Fig. 10 for the same even iteration n = 60. 
Despite the rather different level spectrum for large n in 
the flow diagram in Fig. 9, the actual spin-resolved fixed 
point spectrum is qualitatively very similar to the non- 



magnetic case in Fig. 8. Aside an overall tilt of the level 
structure, all degeneracies and level positions of the lower 
part of the energy spectrum in panel (a) are again fully 
described by elementary single-particle excitations. The 
underlying reason for this similarity of the fixed points 
spectra in the magnetic and non-magnetic case is that, 
apart from the (screened) impurity spin, the system is 
well described by an effective Fermi-liquid picture. With 
the low energy fixed point spectra well reflected in the 
entanglement spectra, a similar tilt in the level structure 
is also observed in Fig. 10(b) when compared to Fig. 8(b). 
Note, for example, that to the lower left of the spectrum 
the same [1-2-1], as well as the [2-4-2, 2-4-2] state se- 
quence with increasing energy is seen. 



IV. SUMMARY AND OUTLOOK 

The reduced density matrices of the NRG by tracing 
out the low-energy sector have been analyzed in detail. 
The low end of their eigenspectra was used to estimate 
the discarded weight gP. in Eqs. (19-21) as a quantita- 
tive and site-resolved measure of the accuracy within the 
NRG. While, in principle, the same reduced density ma- 
trices could also be utilized as the basis for an altered 
truncation criteria similar to the DMRG, this, however, 
requires sufficiently large Mk to start with. In practice, 
this is sufficiently close to a truncation with respect to an 
energy cutoff Ek- Either way, all of this can be easily and 
quickly checked using the proposed analysis in terms of 
the discarded weight which provides a useful quantitative 
tool. 

Furthermore, the dominant correlations of the reduced 
density matrices were analyzed in terms of their entangle- 
ment spectra. Due to the NRG flow towards small energy 
scales, these spectra can be combined into entanglement 
flow diagrams. There different physical regimes can be 
identified similar to the standard NRG energy flow di- 
agrams. Considering that the entanglement spectra are 
obtained solely based on the wave function, the agree- 
ment of the low-energy fixed point spectra are stunning. 
A possible larger disagreement at higher energies, i.e. for 
earlier Wilson shells, on the other hand, depends on the 
specific physical situation. Given the NRG background, 
as an outlook this appears to suggest the following. For 
all energy shells (iterations) n where the entanglement 
spectrum is quantitatively comparable to the NRG en- 
ergy spectrum for the lowest set of states, the reduced 
density matrices themselves are not crucially important 
in the description of the system. Instead, they may be re- 
placed by thermal density matrices in the NRG eigenba- 
sis. In a sense, by tracing out the low-energy sector, the 
resulting reduced density matrices maintain an approx- 
imate thermal character, with implications to thermal- 
ization at a given energy shell. 23 For energy shells with 
a qualitative difference between the energy and entan- 
glement spectra, however, the reduced density matrices 
are crucially important to capture the correct physics in 
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the NRG calculation that explicitly uses data from such 
energy shells. 

A detailed analysis of the deeper connection and the 
explicit differences between the energy and the entangle- 
ment spectra appears interesting, yet is out of the scope 
of this paper. In particular, it also appears instructive 
to analyze the entanglement spectra for non-fermi liq- 
uid systems such as the symmetric two-channel Kondo 
model, as the analysis presented in this paper suggests 
a strong physical connection of the entanglement spectra 
to the underlying physics. 
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